THE SIZE OF RANDOM FRAGMENTATION TREES 



SVANTE JANSON AND RALPH NEININGER^ 

Abstract. We study a random fragmentation process and its associated ran- 
dom tree. The process has earher been studied by Dean and Majumdar |7|, who 
found a phase transition: the number of fragmentations is asymptotically normal 
in some cases but not in others, depending on the position of roots of a cer- 
tain characteristic equation. This parallels the behaviour of discrete analogues 
with various random trees that have been studied in computer science. We give 
rigorous proofs of this phase transition, and add further details. 

The proof uses the contraction method. We extend some previous results for 
recursive sequences of random variables to families of random variables with a 
continuous parameter; we believe that this extension has independent interest. 



1. The problem and result 

Consider the following fragmentation process [7j . Fix b >2 and a random vector 
V = {Vi,...,Vh). Note that the definitions and results below depend only on 
the distribution of {Vi, . . . ,Vb), so it would be more precise to say that we fix a 
distribution on M^; we find it, however, more convenient to state the results in 
terms of a random vector. We assume throughout the paper that < Ij- < 1, 
j = l,...,b, and 

E^^ = i' (1-1) 

i.e., that {Vi, . . . , Vb) belongs to the standard simplex. For simplicity we also assume 
that each Vj < 1 a.s. We allow Vj = 0, but note that, a.s., < V^- < 1 for at least 
one j. 

Starting with an object of size x > 1, we break it into b pieces with sizes 
Vix, . . . , VbX. Continue recursively with each piece of size > 1, using new (indepen- 
dent) copies of the random vector (Vi, . . . , Vb) each time. The process terminates 
a.s. after a finite number of steps, leaving a finite set of fragments of sizes < 1. We 
let N{x) be the random number of fragmentation events, i.e., the number of pieces 
of size > 1 that appear during the process; further, let Ne{x) be the final number 
of fragments, i.e., the number of pieces of size < 1 that appear. 
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This model has been studied by Dean and Majumdar 7 , who found (without 
giving a rigorous proof) that the asymptotic behaviour of N{x) as x — > oo depends 
on the position of the roots of a certain characteristic equation; we give a precise 
version of this in Theorem II .31 below. Some special cases have earlier been studied 
by other authors, see Sectional 

It is natural to consider the fragmentation process as a tree, with the root rep- 
resenting the original object, its children the result of the first fragmentation, and 
so on. It is then convenient to let the fragmentation go on for ever, although we 
ignore what happens to pieces smaller than 1. Let us label each node with the size 
of the corresponding object. 

We thus consider the infinite rooted 6-ary tree Tf,, whose nodes are the strings 
J = ji' ' ' jk with ji G {1, •••,&} and /c > 0. Let B* denote the set of all such 
strings, and let {V^'^\ . . . , V^'^''), J G be independent copies of V. Then node 
J = ji " '3k gets the label xj^^^^^ V^^^ is-i)^ Thus N{x) is the number of nodes 
with labels > 1, i.e. 

By the recursive construction of the fragmentation process, we have N{x) = 
for < X < 1 and 

b 

N{x) = l + ^N^^\Vjx), x>l, (1.3) 
where N^^\-) are copies of the process A^(-), independent of each other and of 

Remark 1.1. Let 'T{x) be the subtree of T{co) = Tf, consisting of all nodes with 
labels > 1. Then A^(x) = |T(x)|, the number of nodes in T{x). More precisely, 
we call these nodes internal nodes of T{x), and we say that a node in T(oo) is an 
external node of T{x) if it has label < 1 but its parent is an internal node. 

Thus A^(x) is the number of internal nodes, and iVe(x) is the number of external 
nodes. Since each internal node has b internal or external children, we have, for 
X > 1, N{x) + Ne{x) = l + bN{x), OT Ne{x) = (6-l)7V(x) + l. Hence the results for 
A^(x) immediately yield similar results for Ne{x) and the total number of external 
and internal nodes A^(x) + iVe(x) too. 

In this paper we thus study the size of the fragmentation tree T{x). Of course, 
it is interesting to study other properties too, such as height, pathlength, .... 

Note that we may define T(x) for all x > simultaneously, using the same V^'''^^ 
this defines {T{x))x>o as an increasing stochastic process of trees. (Equivalently, 
we can label all nodes of T(oo) as above, starting with 1 at the root, and then keep 
all nodes with labels > 1/x). 
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Remark 1.2. We assume for convenience that each object is spHt into the same 
number b of parts. Our method applies also to the case of a random number of 
parts. Indeed, if the number of parts is bounded, we can use the results below 
with b large enough, setting the non-existing Vj := 0. If the number of parts is 
unbounded, we can, under suitable assumptions, use the proofs below with minor 
modifications. We leave this extension to the reader. 

Our main result is Theorem M.IM below on the asymptotic distribution of N(x), 
together with the corresponding estimates for mean and variance given in Theo- 
rem 

We define (with 0^ := 0), at least for Rez > 0, 

b 

c^iz):=Y,^Vf, (1.4) 

and note that (j){z) is bounded and analytic in the open right half-plane {z : Kez > 
0}. Since we assume clearly (/)(1) = 1. Since further < V,- < 1 a.s., the 

function (j)(z) is decreasing for real z > 0; hence (j){z) > 1 when < z < 1 and 
(t){z) < 1 for 1 < 2 < oo. Further, \4>{z)\ < Y.j^\^f\ = (/"(Rez), so \(t>{z)\ < 1 when 
Rez > 1. 

A crucial role is played by the solutions to the characteristic equation 

m = 1. (1.5) 

By the comments above, A = 1 is one root, and Re A < 1 for every root A; further- 
more, there is no real root in (0, 1). 
We further define 

b 

a :=-(!)' {l)=Y,n-VjlnVj), (1.6) 
i=i 

the expected entropy of {Vi, . . . , 14). 

We need a (weak) regularity condition on the distribution of (Vi,...,Vfe). We 
find the following convenient, although it can be weakened to Condition B((5) in 
Section [21 for suitable 5. For examples where this regularity and Theorem 11.31 fail, 
see Example 18.11 

Condition A. Each Vj has a distribution that is absolutely continuous on (0,1), 
although a point mass at is allowed. 

Note that there is no condition on the joint distribution. In one case, however, 
we need a condition including the joint distribution too. (Note that both conditions 
are satisfied if V has a density on the standard simplex, i.e. if (^i, . . . , Vb_i) has a 
density.) 
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Condition A'. The support of the distribution of V on the standard simplex has 
an interior point. 

If Condition A holds, then, by Lemmas 12.21 and 12.11 below, there is only a finite 
number of roots of 0(A) = 1 in {A : Re A > 5} for any 6 > 0. We may thus order 
the roots with ReA > as Ai,A2,...,Am with decreasing real parts: Ai = 1 > 
Re A2 > Re A3 > . . . ; we will assume this in the sequel. If Ai = 1 is the only root 
with ReA > 0, we set A2 = —00 for convenience. 

We let Ai"^ denote the space of probability measures on C, and let 

M^(,):={,,M- : /|.P</,(.)<oo, a„d =7}, 7 e C. 

We let T denote the map (assuming A2 7^ — cxd) 

T:M^^M^, r/^£ , (1.7) 

where (Vi, . . . , V;,), Z^^\ . . . , Z^^^ are independent and C{Z^'^^) = r/ for r = 1, . . . , 6. 
Note that T maps (7) into itself for each 7, since A2 satisfies (/>(A2) = 1. 

We state our main result. The constant a > is defined in (|1.6|) above and f3 is 
given explicitly in Theorem 13. 11 The £2 distance between distributions is defined in 
Section |3 

Theorem 1.3. Suppose that Condition A holds. Then we have: 

(i) // ReA2 < 1/2 then EiV(x) = a'^x + o(^), Var A^(x) ~ Px with /3 > 
and 

^^"^-""'" ^AA(0,/3). 
^x 

(ii) // ReA2 = 1/2 and each root Aj with ReAj = 1/2 is a simple root of 
0(A) = 1, and further Condition A' too holds, then¥,N{x) = a^^x+0{^/x), 
Var(N{x)) ~ Pxlnx with /? > and 

N{x) — a~^x d 



Va^ln: 



^f{o,^3). 



(iii) // Re A2 > 1 /2, and A2 and A3 = A2 are the only roots of 1)1.51) with this real 
part, and these roots are simple, then KN{x) = a^^x + Re(7x'*'^) + 0{x'^), 
with > 0, 7 G C \ {0}, 1/2 < k < Re A2 and 

^^^^^=^,Re(Heii-^^"^-)) = O (^x^~^'^^^ 

for some complex random variable H. Furthermore, >C(H) is the unique fixed 
point of T in A42(7). 
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Remark 1.4. We can regard our process as a general (age-dependent) branching 
process |16| Chapter 6], provided we make a logarithmic change of time as in Sec- 
tion 13 Indeed, there are two versions. For internal nodes, the individuals in the 
branching process live for ever, and give birth at times — In Vi, . . . , — In V^. For 
external nodes, we have a splitting process where each individual when it dies gives 
birth to new particles with life lengths — In Vi, . . . , — In V^. For both versions, we 
obtain a super-critical branching process with Malthusian parameter 1, but the 
identity (jLlj) causes the asymptotics for moments and distributions to be quite 
different from typical super-critical branching processes. 

Remark 1.5. If A2 in (|1.7jl is real, then the stable distributions of index I/A2 are 
fixed points of T. Note, however, that in our case, A2 never is real. Moreover, the 
fixed points we are interested in have finite variance, and are thus quite different 
from stable distributions. 

For the related Quicksort fixed point equation. Fill and Janson found a 
complete characterization of the set of fixed points; in that case, all fixed points are 
formed by combining certain stable distributions with the unique fixed point with 
mean and finite variance. 

Remark 1.6. Condition A' is needed only in part (ii), and is needed only to exclude 
the possibility that for each root Aj with ReAj = 1/2, there is a complex constant 
Ci such that 



This is easily seen to be impossible if Condition A' holds, and even otherwise it 
seems highly unlikely for any particular example, but it seems very likely that there 
are examples satisfying Condition A where V is concentrated on a curve, say, such 
that (dHl) holds. 

We will prove the statements on mean and variance, with further refinements, 
in Section |31 To prove convergence in distribution, we will use a continuous time 
version of the contraction method. We develop a general theorem, that we find to 
be of independent interest, in Section 13 This theorem is applied to our problem in 
Sectional Some examples are given in Sections El and |H1 

Acknowledgements. This research was initiated and largely done during confer- 
ences in Oberwolfach, Frankfurt and Vienna in August and September 2004. We 
thank Luc Devroye, Jim Fill, Allan Gut and Hsien-Kuei Hwang for valuable com- 
ments. 
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2. Further preliminaries 

We define (again with 0^ 
i;{z,w):=CoJyW,yvr'\ = E\ywyV,-\ -cp{z)cp{w). (2.1) 




In particular, iIj{z, z) = K\Ylj=i ~ 4'iz)\ ^ 0) with equahty only if X]j=i — 
(j){z) a.s. 

For Rez,Rew > 0, we have < 1 and thus \i{j{z,w)\ < 26^. 

We say that V is lattice if there exists a number r with < r < 1 such that every 
Vj G {?^"}n>o U {0} a.s.; otherwise V is non-lattice. Basic Fourier analysis applied 
to the probability measure u defined in 1)3. 16() shows that V is non-lattice if and 
only if A = 1 is the only root of (|1.5|) with Re A = 1. (Otherwise, there is an infinite 
number of roots with Re A = 1.) We will assume this, and more, below. 

We introduce a family of regularity conditions that are weaker than Condition A. 

Condition B((5). (Here (5 is a real number with 6 > 0.) 

limsup |(^((5 + it)\ < 1. 

t— >oo 

Lemma 2.1. // Condition B((5) holds for some 6 > 0, then Condition B{6') holds 
for every 6' > 6 as well; moreover 

limsup \ (p{z)\ < 1. 

Rez>5 
Im 2— >oo 

Proof. Choose first e > such that limsup^^o^^ \ (j){5 + it)\ < 1 — 2e, and then A 
such that \(p{5 + it)\ < 1 — 2e if t > ^, and thus also if t < —A. Recall further that 
\(j){6 + it)\ <b for all t. Since 0(z) is analytic, and thus harmonic, in the half-plane 
Tis := {z -.Rez > 6} and bounded and continuous in Tis, (j) is given by the Poisson 
integral of its boundary values |13l Lemma 3.4]: 

/oo 
P^^siy-t)(p{S + it)dt, x>6, (2.2) 
-oo 

where Px{y) = x/{tt{x'^ -\- y^)), the Poisson kernel for the right half-plane. Let 
uj{x + iy) := J^^ P-x-siv — t) dt, the harmonic measure of [5 — lA, 5 + lA]; then 1)2. 2(1 
implies 

/oo 
Px-siy - tM{S + it)\ dt < bio{x + iy) + 1 - 2e. (2.3) 
-oo 

It is well-known, and easy to see, that the set B := {z £ Tig : uj{z) > e/b} is 
bounded; in fact, it is the intersection of Tig and a circular disc jl31 p. 13]. Thus, 
Ai := supjlmz : z £ B} < oo, and if Rex > 6 and \y\ > Ai, then uj{z) < e/b and 
(tO) vields \6(x + iv)\<l-E. □ 
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Lemma 2.2. // Condition A holds, then Condition B((5) holds for every 5 >0. 

Proof. We have E Vj^ = E(e'*^'^^^ l{v^>o}) , the Fourier transform of the distribution 
of InVj- (ignoring any point mass at 0), so by Condition A and the Riemann- 
Lebesgue lemma, EVj^ ^ as t — > oo for every j, and thus (p{it) — > as t ^ oo. 
Hence, Condition B(0) holds, and the result follows by Lemma l2.ll □ 

Lemma 2.3. If Condition B((5) holds for some 5 > 0, then there is only a finite 
number of roots to (/'(A) = 1 with Re X> 6. 

Proof. By Lemma 12.11 all such roots satisfy | Im A| < C for some C < oo. Fur- 
thermore, all roots satisfy Re A < 1, so if further Re A > A belongs to a compact 
rectangle K in the open right half-plane. Since, (p{z) — 1 is analytic and non-constant 
in this half-plane, it has only a finite number of roots in K. □ 

In particular, by the comments above. Condition B{5) with 5 < 1 implies that V 
is non-lattice. 

3. Mean and variance 
We let A denote the set of solutions to the characteristic equation (|1.5j) . i.e. 

A := {A : 0(A) = 1}; (3.1) 

we further define its subsets 

A(s) := {z G A : Re(z) = s}. (3.2) 

In general, 0(A) is defined only for ReA > 0, and we consider only such A in (|3.1() . 
However, in cases where 4> extends to a meromorphic function in a larger domain 
(for example, when is rational), we may include such A too in A; this makes no 
difference in Theorem 13.11 (In Theorem 13. 4| we include all roots in the complex 
plane.) We will use A(s) only for s > 0, where there is no ambiguity. 

Let m{x) := KN{x) and cr^(x) := Var A^(x). We will show the following asymp- 
totics. 

Theorem 3.1. Assume that Condition B((5) holds with < 5 < 1, and let 
Ai, . . . , Xm be the elements of {X £ A : Re A > 5}, ordered with Ai = 1 > Re A2 > 
Re A3 > • • • . Then, the following hold as x ^ 00: 

(i) m(x) ~ a~^x. 

(ii) If further (f)'{Xi) 7^ for i = 1,...,M, i.e., each Aj is a simple root of 
0(A) = 1, then, more precisely, for every 6' > 6, 

M 
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(iii) If 5 < 1/2 and either M = 1 or Re A2 < 1/2, then cr^(x) ~ Px, with 

/9 = a" — / , , ,01 ^-f^ '-^du e (0,00). (3.4 

27ry_^ |l/2 + m|2|l- 0(1/2 + m)|2 ^' ^ ^ ^ 

(iv) If M > 2 and ReA2 = 1/2, and each Aj with ReAj = 1/2 is a simple root 
of 4>{\) = 1, then (T^(x) = Pxlnx + o(xlnx), wii/i 

AeA(l/2) I ^ ^ ^1 

7/, moreover, Condition A' holds (or, more generally, for some Aj G A(l/2), 
(|1.8j) does noi hold), then (3 > 0. 

(v) If M >2 and ReA2 > 1/2, and each Aj with Re A, = ReA2 is a simple root 
of (p{X) = 1, then 

<y\x)= y , . , -, ^^(Ai,Afe)x^»+^^+o(x2i^^^=^). 

Remark 3.2. It follows from the proof that for (i) we do not need Condition B((^); 
it is enough that V is non-lattice. 

Remark 3.3. The case when some (/>'(Aj) = is similar; now terms x^* Inx (and 
possibly x'^* In'^x, d>2) will appear in (|3.3|) . We leave the details to the reader. 

If (/> is a rational function, then ()3.3() can be improved to an exact formula. 
Furthermore, in case (iii) of Theorem 13.11 we then can give an alternative formula 
for /3. 

Theorem 3.4. Assume that (j) is (i.e., extends to) a rational fumction, and let 
Ai, . . . , Ajvf he the roots of (j){X) = 1 in the complex plane. Suppose further that all 
these roots are simple. 



(i) Then 



M ^ ^ 
m(x) = V x^' - — — -, X > 1. (3.6) 



^-A,</.'(A,) 0(0) -1' 



(ii) Assume further that Ai = 1 and ReAj < 1/2 for i = 2,...,M, and that 
Vj > a.s. for every j. Define, for notational convenience, Xq := 0, oq := 
-1/(6- 1) and Ui := -l/(Ai0'(Ai)) fori = 1,...,M. Then a'^{x) ~ /?x. 
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with 

^ = E i-T-X {iZ^^l'^i'"^^^ ^ Vi)'-'^-'" - 20(1 - Afc) + 1 

i,k^l * j,l=l 

M b 
di 



- E >r( E ny-'v^-^' - vi)i{v,<v,} - 0(1 - A.) + 1 

i=2 * j,l=l 

b 

- 2a-^ao ( J] E l^(ln V, - In Vi)l{v^<v,} " «) 

b 

+ a-^(^^E{Vj AVi) -1^ -a-\ 



The proof of these theorems wiU occupy the remainder of this section. We first 
show that all moments of N{x) are finite. 

Lemma 3.5. For every m > 1 and x > 0, EA^(x)™' < oo. Furthermore, 
suPo<y<:,EA^(y)" < oo. 

Proof. For a string J = ji---jk £ S* we denote by |J| = k the depth of the 
corresponding node in T^. Note that we have \{J & B* : < |J| < < b^~^^. 
Hence, if N{x) > b^~^^ for some k > then by (|1.2|) there exists a J = ji • • • j'^ ^ B* 
with vj^^ "''"^^ > 1. Markov's inequality implies that for all > 1 



¥{N{x) > 6^+^) < 



^J€B*:\J\=k U=l ) / 



JeB*:|J|=fc \j=i / 

JeB*:|J|=fc i=l 
= xV(g)^ 

Hence, for all y > 6 we obtain with k = [log^ yj — 1 and 0(g) < 1 that 

nNix) >y)< x^iqf < xV(g)'°*^'^^"^ = y^°SbHg)_ (3.7) 

nqr 

We have (p{q) ^ as g — > 00 since V^- < 1 a.s. and by dominated convergence. 
Hence, for all m > 1 there exists a q > with log^ 0(g) < —m. The tail bound 1)3. 7() 
thus implies EA^(a;)'" < 00 for all m > 1 and all x > 0. 

The final statement follows because < N(y) < N{x) when < y < x. □ 
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We find it convenient to switch from multiplicative to additive notion. We there- 
fore define 

:= -InFj G (0,oo], j = l,...,b, 
N^{t) := N{e^), -oo<t<oo. 

The definition H1.2|) and the recursive equation H1.3() thus translate to 



JeB* 



N^{t) = l + Y^ Ni^\t - Xj), t > 0, (3.9) 

where n!;''\-) are independent copies of the process A^*(-), and N^:{t) = for — oo < 
t < 0. Further define 

:=EiV,(t) = m(e*), 
al{t) := VarA^,(t) = (j2(e*). 
Thus m^:{t) = C7^(t) = for t < 0. Taking expectations in 1)3. 9|) we find 

b 

m*(t) = 1+E^m*(t-Xj), t>0. (3.10) 
i=i 

Let be the distribution of Xj on (0, oo); this is a measure of mass 1 — P(Vj' = 0); 
let further fi := Yl^j=i f^j- Then 1)3. 10() can be written 

b 

m^{t) = l + '^lj,j*m^{t) = 1 + fi*m^{t), t>0, (3-11) 
i=i 

where fj, * f{t) = f{t — x) d^{x). This is the standard renewal equation, except 
that /X is not a probability measure. 

Similarly, conditioning on Xi, . . . , X^^ for t > 0, 

E((iV,(t)-m,(t))2 I Xi,...,Xf,) 

b 



2 

, . . . , Xfe 



E ( Ni^^ {t-Xj) + l-m^{t) 

6 

Var Ni^\t - Xj) + l-m^{t)\Xi,..., X^ 
i=i 

6 , 
+ (^^m=,(t-Xj) + l-m*(t)y 
j=i 

^ a2(t - Xj) + m*(t -Xj)- m,{t) + 1 
i=i i=i 
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Taking the expectation we obtain 

b 

alit) = EY,(^l{t-Xj) + h{t) = fi*al{t) + hit), t>0, (3.12) 



where, recalhng (|3.1U|) . 

b , 
h{t) := E(^^m*(t - Xj) - m^t) + ij 



= E(^m,(t-Xj) -m,(t)) +2(E^m,(f-Xj) -m,(t)) +1 

i=i j=i 
b 2 
= E(j2m,{t- Xj) -m4t)^ -1. (3.13) 

i=i 

Both ()3.1ip and (|3.12p are instances of the general renewal equation (|3.14p below, 
and from renewal theory we get the following result. We say that a function on [0, oo) 
is locally bounded if it is bounded on every finite interval. 

Lemma 3.6. Assume that V is non-lattice. Let f be a locally bounded measurable 
function on [0, cxd). Then the renewal equation 

F = f + n*F (3.14) 

has a unique locally bounded solution F on [0, co). If further M is non-lattice, then 
we have the following asymptotical results, as t ^ oo, 

(i) // / is a.e. continuous and f*{t)dt < oo, where f*{t) := 
sup„>te-"|/(M)|, then F{t) = (7 + o(l))e*, with j = a'^ f{t)e-^ dt. 

(ii) // f{t) = e*, then F{t) ~ a^heK 

(iii) // f{t) = e^* with Re A = 1 and Im A 7^ 0, then F{t) = o{te^). 

(iv) // f{t) = e^* with ReA > 1, then F{t) ~ (l - 0(A)) "^e^*. 

Proof. For a function / on (0, 00) and z G C, we define, when the integral exists, 
the Laplace transform f(z) := e"^^ f{t) dt. Similarly, the Laplace transform of 
fi is 

POO b ^ 

Ji{z) := / e-'^df^{t) = ^Ee-^^^- = ^ET/^ = <p(z), (3.15) 

i=i j=i 

at least for Re z > 0. 

Since fi is not a probability measure, we define another ("conjugate" or "tilted") 
measure u on [0, 00) by 

du{x) = e-"" dn{x). (3.16) 
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Then is a probability measure because, by 

poo ^ poo ^ ^ 

i/[0,oo) = / e-^d^(x) = V / e-^(i/ij(x) = VEe"^^ = VeT/^ = 1. 
Further, the mean of the distribution v is 

<>oo j-OO ^ ^ 

Kiy= xdv{x) = / xe-^(i^(x) = VE(Xje-^^) = ^^{{-\nVj)Vj) = a 
Jo Jo ^.=1 

(3.17) 

and the Laplace transform is, for Kez > 0, recalling f|3.15p . 

/"CO roo 

Hz):= e-^''dv{x)= e'''-'''d^i{x) = Jl{z + l) = (t>{z + l). (3.18) 
Jo Jo 

Let fi((t) := e"*/(i) and G(t) := e~*F{t). Then translates to 

/>oo 

G{t) = e-*F{t) = e-^f{t) + / e-*F(t - x) d^l{x) 

Jo 

/•oo 

= g{t)+ / G{t- x)e~^ dn{x) = g(t) + u*G{t). 
Jo 

In other words, G satisfies the renewal equation for the probability measure v, so 
we can use standard results from renewal theory. 

First, it is well known that the equation G = g + v * G has a unique locally 
bounded solution which is given by G = Yl'^=o * 5'' t\ms F = Yl^=o /^*" * /> 
see e.g. Theorem IV. 2. 4] (which also applies directly to F). If we let Yi, Y2, ■ ■ ■ 
be i.i.d. random variables with the distribution v, and let Sn '■= ^1^1 this can be 
written 

00 

G{t) = ^K{g{t - Sn)lis^<t}) =Ej29{t- Sn). (3.19) 

n=0 Sn<t 

Under the assumptions of (i), /* is non- increasing and integrable; further, 
sup/* < sup[o^i] I/I + /*(1) < 00, so /* is bounded too. Hence [T| Proposition 
IV.4.1(v),(iv)] shows that /* and g are directly Riemann integrable. The key renewal 
theorem ^ Theorem IV. 4. 3] and (|3.17|) now yield G{t) a^^ g{x) dx = 7, 
which proves (i). 

In case (ii) we have g{t) = 1, and thus G{t) ~ a~^t by the elementary renewal 
theorem ^ IV. (1.5) and Theorem 2.4]. 

For (iii), g{t) = e^^~^^^ = e*''* for some real 6 7^ 0. The solution to (|3.14|) may 
be written ^ Theorem IV. 2.4] G{t) = Jq g{t — x)dU{x), where U is the locally 
bounded solution to U = 1 + 1^*11 (i.e., U = G for case (ii)). Since, in analogy with 
p.l7() . / x"^ dv{x) = E((ln V^)^Vj) < co, the distribution v has finite variance, 
and the renewal theorem has the sharper version 01 Theorem XI. 3.1] 

U{t) = a~^t + c + R{t), 
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where c is a certain constant (J dv /2a^) and R{t) ^ as i ^ oo. Hence, using 
integration by parts for one term, 

G{t) = f e'^(*-^)a-^ dx + ce'** + f e"'(*-^) dR{x) 
Jo Jo 

= 0(1) + R{t) - i?(0)e'** + ib f e'^(*-^)i2(x) dx = o{t). 

Jo 

For (iv), we have by (|3.19|) with g{t) = e^^~^^^, using dominated convergence and 

(Em, 



oo 



n=0 

oo oo oo 

n=0 n=0 n=0 

oo 

= j]0(Ar = (i-0(A))-\ 

n=0 □ 

Proof of Theorem VJ.l\ We first apply Lemma I3.6f i') to ()3.11|) . with f{t) = 1 for 
t > 0, and obtain 7 = and m^:{t) ~ a~^e*, which proves Theorem I3.1f i). 

To obtain more refined asymptotics, we use Laplace transforms. Let H{t) := 
l{t>o} (^^6 Heaviside function), and note that H(z) = e~^^ dt = 1/z, Rez > 0. 
Since the Laplace transform converts convolutions to products, the renewal equation 
(|3.11() yields rn^{z) = H{z) + Jl{z)m'^:{z), and thus 

~. X H(z) 1 

l-^J.{z) z{i-<t){z)) 

for z such that the transforms exist. By the estimate m^it) ~ a^^e* above, m^:{t) = 
0(e*) and thus rru{z) exists for Rez > 1. Consequently, ()3.2Up holds for Rez > 1, 
and can be used to extend frh,{z) to a meromorphic function for Rez > 0. 

We want to invert the Laplace transform in (|3.2U|) . This is simple if is rational, 
yielding ()3.6p . (Note that (f)(0) = E|{j : Vj > 0}| > 1.) In general, there are 
difficulties to doing this directly, because rn^,{z) is not integrable along a vertical 
line Rez = s; it decreases too slowly as | Imz| 00. We therefore regularize. Let 
e > 0, and let := H * e~^l[o,e]; thus 

'0, t<0, 
He{t) = ll-t/e, 0<t<e, 
1, t> e. 

Let m^<e = J2'^=o ^^*^ * ^£ locally bounded solution to m^K^ = + fi * m^£. 

Note that Heit) < H{t) < Heit + e), and thus 

m^e^t) < m^{t) < m^eit + e). (3.21) 
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We have 



Heiz)=H{z)e-^ I e-'Ut = - — — , Re2;>0, 



and we find, arguing as for (|3.2U|) above, 

He{z) 



first for Rez > 1, and then for Rez > 0, extending m^^ to a meromorphic function 
in this domain. This function decreases (using Condition B((5) and Lemma l2.1|) as 
|Imz|~^ on vertical lines Rez = s > 5, and is thus integrable there. Hence, the 
Laplace inversion formula (a Fourier inversion) shows that for any s > 1 and t > 0, 



^ ps+ioo 

m*e{t) = — e''f^e{z)dz. (3.22) 



We may, increasing 6 a little if necessary, assume that (/)(z) = 1 has no roots with 
Rez = 5; in cases (iii), (iv) and (v) we may similarly assume that each A G A 
with Re A > 5 has 0'(A) ^ 0. It is then easy to show, using Condition B(5) and 
Lemma 12.11 that we may shift the line of integration in (|3.22|) to Re z = 6 and 
obtain, for < e < 1, 



"/ S — icxD ^ 



^ -Xi(t)'{Xi) eXi 



-t 



Now choose e := e ; we then obtain 



m ' " 



«(') = Erif;^ + o(i) + o(e«(i + *)), t>o. 

Replacing t by t + e, we obtain the same estimate for m^^{t + e), and thus (|3.21|) 
yields 

».(')= Er^ + 0(."'), t>0, (3.23) 

which yields Theorem 13.1( 11). 

For the estimates of the variance, we use Lemma 13.61 and ()3.12|) . It is easily seen 
(by dominated convergence) that h in (|3.13|) is a.e. continuous. Choose 5' > 5 with 
5 <5' < Re Am; in case (iii) with M = 1, let further 5' < 1/2. 
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Note that then (|3.23)) trivially holds for t < too. Hence, 

b M yb^ g(t-x,)x, _ ^tx, 
Y^m^it- X,) - m^it) = Y: + 0{e''') 

j — 1 i — 1 

^ tX ( b \ 

where we use the fact that Ai = 1 and thus Yl^j=i ~ 1 = Sj=i ^3 ~ ^ ~ ^■ 
Consequently, by ^^TT^ and (EH}, letting a2 := ReA2 > 5' if M > 2, and (72 := 6' 
if M = 1, 

''W-gE ,.,;,,(;.)^(,,) ^(A..A.) + 0(.'<^'^-.»)^ (3.24) 

For Theorem EIHiii), yields h{t) = 0{e'^''^^) with c72 < 1/2, and 

Lemma [3.6f i) applies to (|3.12p . yielding cT^(t) ~ 7e*. We postpone the calcula- 
tion of /3 = 7, verifying (|3.4p . to Lemma IT71 

In Theorem l3 ■ 1 T iv) and (v), we treat the terms in (|3.24() separately, using linearity; 
for the error term we also use monotonicity and comparison with the case f{t) = 
gt(<5'+o-2)_ fpj^g results follow by applying Lemma l3.6f iVfiv) to 1)3. 12() . letting t := 
Inx, with the leading terms coming from the case ReAfc = ReAj = ReA2, and for 
(iv) further A^ = Aj. 

Furthermore, by ()3.5() . /3 = in (iv) only if for every Aj G A(l/2), we have 
V'(Ai, Ai) = E I Ej - 'A(Ai)P = 0, and thus (HH) holds. □ 

Lemma 3.7. // ReA2 < 1/2 or M = 1, then, with h{t) as in (fXT^ . 

1 /"°° ^(1/2 + 1^,1/2-1^) 
2^y_^ |l/2 + in|2|l-,/.(l/2 + in)|2 

Proof. Write f{t) := (m,(t) - a"^e*)e"*/2^ -oo < t < oo. Thus, by ^2E^, 
f{t) = 0(e-(i/2-<T2)t) for t > and f{t) = -a'^e^l'^ = 0(e-l*l/2) for t < 0. 
In particular, / E L2(— oo,oo). Furthermore, the (two-sided) Laplace transform 
f(z) := /(^e^*^ is analytic for -(1/2 - 0-2) < Rez < 1/2. 

Define further fi{t) := /(t)e*/2 = m,(t) - a'^e* and /2(t) := /i(t)l{t>o}. Then 
f2{t) = 0(6°^^*), and thus the Laplace transform f2{z) is analytic for Rez > CJ2. 
For Rez > 1 we have, by 

72(2)= / e-^''{m^{t)-a''^e^)dt = m4z)-a~\z-iy^ = — — -; 

Jo z{l-(j){z) a{z-l) 

by analytic continuation, this formula holds for Rez > a2- Consequently, for CJ2 < 
Rez < 1, 

h{z) = h{z)+f e~'^{-a~'e')dt = f2{z)-a-Hl-zr' ^ 



z{l-<P{z) 
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Since f{z) = fi{z + 1/2), we find the Fourier transform 



1 



/(n):=/ e-'^'f{t)dt = f{m) = f\i^ + m) = -^— 

Jo (2 +m)(l-(/'(2 +1^^)) 



. (3.25) 



Next, since e'^' = Ej = 1> 

b b b 

Y^mS- Xj) - mS) = 5^ - Xj) - h{t) + ^ e*" 
j=i j=i j,=i 

b 

= Y^f,(t-x,)-Mt), 

so by p.l3|) . and defining ^'(x, y) := /(t - a;)/(t - y) dt, 

/"oo /"oo poo ^ 

I h{t)e-^dt + l= I {h{t) + l)e-^dt = E (S^ fi{t - Xj) - fi{t)f e'^ dt 

Jo Jo J-OD • ^ 



/oo I ^ 
\j2e-''^/'f{t-X,)-f{t) 



dt 



E ^ e-^^/2"^'=/2^'(Xj-,Xfc) - 2E^e-^^/^^(Xj,0) + ^'(0,0). 



i,fc=i 

By Parseval's relation and f = f, 

^ />oo 



^'(2;,y) 



27r 



e-™^/(u)e-i«y/(u) dn— / |/(n)|2e'"(f-^) dt/. 



27r 



Hence, 



/ h{t)e-'dt + l=E- \f{u)nYe- 
Jo J -00 S-,fc=i 



■Xj/2-Xk/2+iu(Xk-Xj) 



1 /"°° I 



6 



fc=l ^ 



1/2-m 



1 



(in 



E — / |/(n)|2(V'(l/2 + in,l/2-in) + |</.(l/2 + in)-l|2)du. 



Using ^n^ . 

IE—/ |/(n)p|0(l/2 + in)-l|2dn 
and the formula follows by (|3.25p . 



du 



|l/2 + in|2 J_^l + u 



du 



1 -i- „2 



2vr, 
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Since ip{z,z) > 0, and by dominated convergence is continuous for Rez > 0, 
it follows from (|3.4|) that /? = only if for every z with Rez = 1/2, = 
IE I Ej ^/ - </'(-2)P = 0> and thus 

^V,i/2+™ = ^(l/2 + in) (3.26) 
i=i 

a.s., for every real u. Considering first rational u, we see that a.s. (|3.26|) holds for 
all real u. 

However, for any realization (Vi, . . . , Vh) and e > 0, it is by Dirichlet's theorem 
jl4| Theorem 201] possible to find arbitrarily large u with Re > (1 — £)V^^'^ 

for j = 1, . . . , 6, and thus (|3.26j) implies that limsup^^oo \4>{^ + i'")! ^ </'(l/2) > 1- 
This contradicts Condition B((5) and Lemma l2.ll Hence /3 > 0. □ 



Proof of Theorem Yj.4\ (i): As remarked above, 1)3. 6(1 follows by inverting the 
Laplace transform in 1)3. 2U() . using a partial fraction expansion. 

(ii): Note first that </)(z) — > as z ^ +oo (by dominated convergence), and thus, 
(f) being rational, </>(oo) = and (j){z) — > as \z\ — > oo; hence Condition B((^) holds 
for every 6. We thus see that the conditions of Theorem 13.1( 111') are satisfied, and 
from the proof above we see that, with h given by 1)3. 13() . 

/"OO /'OO 

P = a'^ J h{t)e'^ dt = a'^ J h{ln x)x~^ dx. 
We have, by (|3.6|1 . m{x) = Ei^o '^«^'^'-'-{a;>i} and thus 

M 

m^,(lnx — Xj) = m{xe~^^) = in{xVj) =''^^ai{xVj)^'l^^yy-iy 

1=0 ~ ^ 

Hence, letting Vq := 1, Eq = —1 and Sj = 1 for j > I, and recalling ()1.1|) . 

6 Mb 

H{x) := ^m*(lnx - Xj) - m*(lnx) = ^aiX^" {^y^"'^{^>y-^} - l{x>i}) 

3=1 i=0 3=1 

b b 

iT^l i=0 i = l 
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By (Trnri) . this leads to 



/oo /"OO roo 

h{\nx)x-^dx + l = J EH{xfx-'^dx = EJ H{xfx' 

^ poo 

V e^e/K^'V^^*^ / x^'+^^-^dx 



^dx 



b 



b b y-i 

15,^1 j=0 «=i ■'^i ^ ^ ' 

+ Kal V FjV; / dx 

and the result follows by straightforward calculations, noting that ai = a^^ . □ 

4. ZOLOTAREV METRIC AND MINIMAL Ls METRIC 

In this section we collect properties of the minimal Ls metric and the Zolotarev 
metric that are used subsequently. 

We denote by A^*^ the space of probability measures on W^. The minimal Lg 
metric £5, s > 0, is defined on the subspace A^f C A^'^ of probability measures with 
finite absolute moment of order s by 

^s{^^,y) :=inf {||X-y||f 1 : X = /i, y = 1/}, n,u£Mi 

where II^H^ := (ElX^^/^ denotes the Ls norm of X. The infimum is taken over 
all random vectors of X, y on a joint probability space with the given marginal 
distributions // and v. (In other words, over all couplings (X, Y) of and v.) We 
will also use the notation £s{X,Y) := £s{C{X), /:{Y)). 

For s > 1 and 7 G M'^, we denote by A^f (7) C A^^ the subspace of probability 
measures with expectation 7. The pairs (A^^,^^)) s > 0, and {Ais{'y),is), s > 1, 
are complete metric spaces and convergence in ig is equivalent to weak convergence 
plus convergence of the absolute moments of order s. 

Random vectors {X,Y) with X = fi, Y = u, and ls{ii,y) = \\X — Y\\f^^ are 
called optimal couplings of (fJ-,!^)- Such optimal couplings exist for all /U, 1/ G A^^. 
These properties can be found in Dall'Aglio [S], Major 23 , Bickel and Freedman 0, 
and Rachev j^. Similar properties hold for probability measures on (because 
we use corresponding notations. 

The Zolotarev metric Cs, s > is defined by 

CsiX,Y):=a/:{X),/:iY)):= sup | E(/(X) - /(y))| (4.1) 

where s = m + a with < a < 1, m = [s] — 1 > is an integer, and 
J's := {/ G C^iR^R) : \\f^^\x) - /(™)(y)|| < ||x - y\n, 
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where C"'(M'^,M) denotes the space of m times continuously differentiable functions 
/ on M'^ and their mth derivative. 

The expression (!^s{X,Y) is finite if X and Y have finite absolute moments of 
order s and all mixed moments of orders 1, ... ,112 of X and Y coincide. 

The metric (s is ideal of order s, i.e., we have for Z independent of {X,Y) and 
any d x d square matrix A 

UX + Z,Y + Z)< Cs{X, Y), Cs{AX, AY) < \\A\\lp UX, Y), 

where ||^||op := sup||„||=;^ ll^^ll denotes the operator norm of the matrix. Conver- 
gence in (^s implies weak convergence. For general reference and properties of Cs we 
refer to Zolotarev [2^1 EOl and Rachev . 



5. General contraction theorems in continuous time 

In this section we extend a general contraction theorem for recursive sequences 
{Yn)n>o of d-dimensional vectors as developed in Neininger and Riischendorf [211 to 
families (lt)t>o of d-dimensional vectors with continuous parameter t G [0, 00). For 
this, we assume that we have 

K 

Y,^^Ar{t)Y^%+h, t>To, (5.1) 
1 

r=l 

where K is a positive integer and tq > 0, and we have {Y^^^)t, . . . , {Y^^^)t, 
{Ai{t), . . . , AK{t),bt, T'^^^)t independent, where = (T^ \ . . . , T^*^) is a vector 
of random indices tJ*'' S [0, t], the Ar{t) are random dx d matrices for r = 1, . . . , K 
and bt is a random d dimensional vector. Finally, in (|5.1j) . we have that for each 

(r) 

t > 0, Yt and Y^ are identically distributed for all r = 1, . . . , K. 

We assume that all Yt as well as Ar{t), bt and T^*) are defined on some proba- 
bility space {il,J^,n), and that they are measurable functions of {t,uj). (This is a 
technicality to ensure that the sum in (|5.1|) is well-defined. Note, however, that the 
joint distribution of Yt for different t is irrelevant.) 

We introduce the normalized random vectors 

Xt:=C;^^\Yt-Mt), t>0, (5.2) 

where Mt G M"^ and Ct is a symmetric, positive definite square matrix. We assume 
that Mt and Ct are measurable functions of t; further restrictions on Mt and Ct will 
be given in Convention C. The recurrence ()5.1|) implies a recurrence for Xt, 

t>ro, (5.3) 

r 
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with independence relations as in (|5.H) and 



(5.4) 



As for the case with integer indexed vectors we estabhsh a transfer theorem of the 
fohowing form: Appropriate convergence of the coefficients A^r^ A^, 6^*) b* 
impUes weak convergence of the quantities Xt to a hmit X. The distribution C{X) 
of X is a fixed point of the hmiting equation obtained from (|5.3jl by letting formally 



oo: 



K 



X = +6*, 



(5.5) 



r=l 



where (A*, . . . , 6*), X^^), . . . , X(^) are independent and = X for r = 
1,. . . ,K. To formalize this we introduce the map T on the space Ai"^ of proba- 
bility measures on R*^ by 



' K 

+6* 

\r=\ 



(5.6) 



where {A\, . . . , A\,}f), are independent and /:(ZW) = r? for r = 

1, . . . , ii'. Then X is a solution of ()5.5j) if and only if C{X) is a fixed point of T. 

We make use of Zolotarev's metric with < s < 3. To ensure finiteness of the 
metric subsequently we make the following assumptions about the scaling imposed 
in (Ol) : 

Convention C. For 1 < s < 3 we assume that Mi = El^. For 2 < s < 3 we assume 
that Cov(yt) is positive definite for all t > ri with a ri > tq and that Ct = Id^ for 
< t < n and = Cov(yi) for t > n. 

This convention implies that Xf is centered for 1 < s < 3 and has Id^ as its 
covariance matrix for 2 < s < 3 and t > ti. (For < s < 1, Convention C is void.) 

Theorem 5.1. Let < s < 3 and let (^)i>o be a process of random vectors satis- 
fying 1)5. 1|) such that \\Yf\\s < oo for every t. Denote by Xt the rescaled quantities 
in ()5.2() . assuming Convention C. Assume that \\A^^'^\\s < oo, ||6^*^||s < oo and 
supq<„<j ||Xu||s < oo for every t > 0, and 

is 



4\...,A^lb(^^ 

K 



{A\,...,A*K,b* 



r=l 



E 



< 1, 



op 



(5.7) 
(5.8) 

(5.9) 
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for every r > and r = 1, . . . , K . Then Xt converges in distribution to a limit X, 
and 

UXuX)^Q, t^oo, (5.10) 
where C{X) is the unique fixed point ofT given in (|5.(i|) subject to \\X\\s < oo and 



EX = forl<s<2, 
EX = 0, Cov(X) = Idd for2<s<3. 



(5.11) 



Proof. This proof is a continuous extension of the proof of Theorem 4.1 in Neininger 
and Riischendorf j2l] for the discrete time case. The existence and uniqueness of 
the fixed point of T subject to (|5.11|) is obtained as follows: For 1 < s < 3 equation 
(tO) implies E = for all t > 0, thus by we obtain E 6* = 0. For 2 < s < 3 
equation (|5.3() implies that for all t > ti 



Idrf = Cov(Xt) 



E 



tr 



■E 



K 



where 6*'" denotes the transpose of a vector or matrix and Ct := Cov(Xt); recall 
that Ct = Id when t > ti. 

By (|5.7j) . (|5.9|) and Holder's inequality this implies 

K 



E[b*{b*f] +E 



.r=l 



Ida. 



Now, Corollary 3.4 in j24j implies existence and uniqueness of the fixed-point 
Since 

K K 

EEii4*^iiop-E^p;ii^p = e<i 

r=l r=l 

there exist G (^,1) and T2 > ri such that for all t > r2 we have 

K 

EEii4*^ii^p<e+<i. 



(5.12) 



(5.13) 



r=l 



Now, we introduce the quantity 



r=l 



i>Ti, (5.14) 



where (4*\ • • • , ^W), X^, . . . , X^, (xf ^), . . . , (X^^^) are independent 



with X^*^) ~ X and xf ~ Xj for r = 1, . . . , K and t > 0. Comparing with (tO) 
we obtain that is centered for 1 < s < 3 and has the covariance matrix Id^ for 
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2 < s < 3 and t > ti. Hence, Cs distances between Xt, Qt and X are finite for all 
t > Ti. The triangle inequality implies 

A(t) := UXuX) < UXuQt) + UQt,X). (5.15) 

As in the proof for the discrete case we obtain QsiQt^X) — > as t — > 0, where we 
use that supo<t<T-2 II -''^t lis < °o. 

The first summand of ()5.15|) requires a continuous analog of the estimate in the 
discrete case. Using the properties of the C,s metric, we obtain, for t > n, 

K 

r=l I i °P 

and, with (|5.15|) . and rt ■= Cs{Qt,X) it follows 

K 



A(t)<E^l 



op 



A(rW) + n. 



(5.17) 



Now, we obtain A(t) — > in two steps, first showing that (A(t))t>o is bounded and 
then, using the bound, that A(t) 0. 
For the first step we introduce 

A*{t) := sup A('u). (5.18) 

T2<n<t 

We have A*{t) < oo for all t > T2, since, for T2 < u < t, we have CsiXu,X) < 
Csi\\X\\l + \\Xu\\s) < C'sdl-'^lls + sup^^<^<t\\Xu\\l) < oo with a constant Cs > 0, 
using |29| Lemma 2]. By definition. A* is monotonically increasing. With R := 
supj>^2 n < CO we obtain for T2 < u < t, from (|5.17|) . (|5.18|) and (|5.13|) . 

K 



r=l I- J 

< (+A*{t) + R. 



op 



A*{u) + R 



Hence, we obtain A*(t) < i+A*{t) + R, thus A*{t) < R/{1 - ^+). This implies 



R 

A*(oo) := supA(t) < - 

t>T2 1 - 



< oo. 



(5.19) 



For the second step we denote L := limsup^^o^^ A(^)- For every e > there exists 
a T3 > r2 such that we have A(t) < L + e for all t > t^. Thus, from 1)5. 17(1 we obtain 

K K 



op 



{L + e)+ rt 



r=l r'=l ^ ' 

and letting t ^ 00 we obtain by (|5.9)) and (|5.12)) 

L < C(L + e). 

If L > 0, this is a contradiction for < e < L[l — Hence, we have L = 0. This 
proves ()5.1U|) . Finally, recall that convergence in Q implies weak convergence. □ 
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As a corollary we formulate a univariate central limit theorem that corresponds 
to Neininger and Riischendorf |24| Corollary 5.2] for the discrete time case. For this 
we assume that there are expansions, as t — > oo, 

¥.Yt = f{t) + o{g^'\t)), \8.r{Yt) = g{t) + o{g{t)) (5.20) 

with functions / : [0, oo) ^ M, (7 : [0, cxd) — > [0, 00), with 

sup \ f{u)\ < 00 for every t > 0, lim g{t) = 00, sup5((n) = 0{g{t)). (5.21) 

u<t u<t 

Thus, for some constant C > 1, g{u) < Cg{t) when < n < t. 
Then the following central limit law holds: 

Corollary 5.2. Let 2 < s < 3 and let Yt, t > 0, be given s-integrable, univariate 
random variables satisfying (|5.1|) with Ar{t) = 1 for all r = 1,. . . ,K and t > 0. 
Assume that sup^<^E|yu|* < 00 for every t, and that the mean and variance of Yt 
satisfy (|^^ with ((??7T|) . //, as t ^ 00, 



9{t) '■■■'V 9it) 



{AI,...,A*k), (5.22) 



^ [bt - fit) + f; /(tW) j A, 0, (5.23) 



and furthermore 



Y.{A;f = i, f(\J{a; = i}\ <i, (5.24) 

r=l \r=l / 



then 



(5.25) 

Proof. We begin by replacing g{t) by max{g{t),l); by 1)5. 21(1 . this does not affect 
g{t) for large t, and it is easy to see that (tPT^ . (TTTTI) . (tPT^ . (tPT^ still hold. We 
may thus assume that g{t) > 1 for every t. 

Denote Mt := EYt and a"^ := Yai{Yt). By (E^Ol), (^t/dit) ^ 1- (Ah unspecified 
limits are as f — > oo.) Choose ti > tq such that \g{t) < (t| < 4g((t) for t > ri. Let, 
as in Convention C, Ct := 1 for t < ti and Ct '■= (t^ for t > ti, and write cJt := C/^'^ 
and£(t) :=at/5(t)i/2_i = (c'tMt))'/'-!. For t > n, e(t) = (Var y^Mt)) - 1, 
so by dOni), 

e{t) ^0 as t ^ oo. (5.26) 
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Further, Ct/9{t) = l/g{t) < 1 for t < n, while Ct/g{t) = a^/g{t) < 4 for t > n. 
Hence |e(t)| < 1 for all t. With ((OJ and Ar{t) = 1 we have, for t > ti, 



(5.27) 
(5.28) 



Since ^(rr^*^) < Cg{t) by dHUIl), we have, for t > n. 



/l(i) 



9{Tn 



W\l/2 



0-t 



e(TW) 



,5(T, 



(t)^l/2 



0-t 



< sup 

u<t 



e(n) 



1/2 



0"t 



(5.29) 



For any 6 > 0, there exists, by 1)5. 26() . r((5) > ri such that |e(t)| < (5 when t > t{5). 
Thus, if t{6) <u<t, then 



e(u) 



9{u 



)l/2 



0-* 



On the other hand, if n < t{5), then 



e(n) 



,1/2 



^ Cg{T{6)Y/^ ^ 

~ CTt 



as t — > oo. Hence, sup„<j \£{u)g{uy^'^ /at\ < 2C6 for sufficiently large t. Since 6 > 
is arbitrary, it follows that the right hand side of (|5.29j) tends to as t — > oo, and 
thus (f5:29j) yields 



9{Tr 



0. 



(5.30) 



Since g{tY/'^/(Jt 1, yields g{Ti^^^f/'^/at A;, which combined with (OHl) 

yields A^f^ A*, jointly for r = 1, . . . , fc. 

Next, for any e > 0, there exists by (fOI)) > n such that |M(-/(t)| < eg{ty/'^ 

t > Tgr. Consequently, if T^*^ > r^, then 

|M - /(tW)| < £5(^^)1/2 < Ceg(f)V2. 



Since sup„<^^ |M„| and sup„<^^ 1/(^)1 finite, the same estimate holds for T^^ < 
too, provided t is large. Consequently, |M (t) — /(Tj*^)|/g(t)^/^ < Ce if t is large 
enough. It follows that || M^(t) - /(Tj*^)||,/5(t)i/2 ^ as t ^ oo, so by 

(t^:^ and 0. 

We apply Theorem 15 . II with 2 < s < 3; we have shown that (|5.7() holds with b* = 
0. The two assumptions in (|5.24|1 and s > 2 ensure that we have lE^^;^ l^rl'' < 1- 
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Finally, by (|5.H()|) . for every r and r, 



< 



< 



9{Tr 



it)\l/2 



+ 



+ o(l) ^ 0. 



Now, Theorem 15.11 implies (Yt — Mt)/(Jt X, where CiX) is characterized by 
\X\\s < oo, EX = 0, Var(X) = 1, and 

K 

X=Y,KX^"\ (5.31) 



r=l 



with assumptions as in H5.5(l . Since 

Er=i(^r)^ = 1 this is solved by C{X) = 1). 
Consequently, 



Yt - Mt ^ 



AA(0,1), 



which, in view of H5.20() . implies the assertion. 



□ 



The following theorem covers cases where the previous central limit theorem of 
Corollarv 15.21 fails due to the appearance of periodic behavior. For this we assume 
that there is an expansion of the mean, as t ^ oo, 

Eyt = /(t)+Re(7t^) + o(r), 



with a function / : [0, oo) 
denote 



(5.32) 



7 G C \ {0}, and A e C with o := Re(A) > 0. We 



it)' 



r = 1 



, . . . , -t i , 



(5.33) 
(5.34) 



Note that Ar^ in general is complex, while 6^*^ is real. 

Theorem 5.3. Let Yt, t >0, be given square-integrable, univariate random vari- 
ables satisfying 1)5. 1|) with Ar{t) = 1 for all r = 1,. . . and i > 0. Assume that 
sup„<(E|yup < oo for every t > and that the mean of Yt satisfies (|5.32() with 
A = cr + ir and a > 0, and some locally bounded function f{t). If, as t ^ oo, 

{Af\. . . , A^) ^ {Al, ...,A*j,) and \\b^'^h - 0, 

and furthermore 



(5.35) 



K 



e^|a;|2<i, 



(5.36) 



r=l 
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then, as t oo, 

£2 ( ^'~J^^\ Re (^Xe'^^'^'j ) ^ 0, (5.37) 
where C{X) is the unique fixed point in (7) of 



c{^A;Z^^)^ , (5.38) 



where {Al, . . . , A*j^) , Z^^\...,Z^^^ are independent and L{Z^^^') = r/ for r = 
1,...,K. 

Proof. We extend an approach based on the contraction method from Fill and Kapur 
We may assume that tq > 1. 
First, for technical convenience we show that we further may assume Yt = and 
f{t) = for < t < 1. Let {Y}'"'^*)t, r = l,...,K,he another set of copies of {Yt)t, 
independent of each other and of everything else. We may replace Y^^^ in 1)5. 1(1 
by y/^^l{(>x} + which has the same distribution and independence 

properties. Hence Y^^^l^^^^ satisfies (jS.ip (for t > tq > 1) with bt replaced by 
bt ■■= bt + Er^^w* l{r(*'<i}- ™s replaces 6^ by feW with 



I J f- 



^{r(*'<i} 



so II^W - 6^112 = Oit-") and (f^:^ still holds. We may thus consider y/''^l{i>i} 

(r) 

instead, and thus we may assume that 1^ = when t < 1. Similarly, we may 
assume that f{t) = for t < 1, changing 6^*-' by 0{t^'^). 
With Xt := {Yt - fit))/^ for t > and Xq := we obtain 

^* = EH^ t>ro, (5.39) 



r=l 



with as given in 1)5.34(1 . 

Next we prove that the restriction of T defined in (|5.38l) to 7^2(7) 
maps into 7^2(7) is Lipschitz in £2 with Lipschitz constant bounded by 

Note that (|5.36l) implies ||^*||2 < 00 for all r = 1, . . . , This implies that T{ri) 
has a finite second moment for all rj G A^2- Next we claim that Er^i^^* — ^■ 
This implies that T(r]) has mean 7 for all r] G ^^2(7)- To prove Er^i^^* ~ 1; 
note that ()5.32|) implies EXt = Re(7t''^) + o(l) as t ^ 00. On the other hand, the 
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right hand side of (l.^.lfflj) has mean, using E6W ^ 0, 



K 



r=l 



Re(7(TW)-) 



+ o{l) = RehY^E^-^ \ +0(1) 
= Re (^rYlEA*}j +o(l 



where we also used that E(rJ*Vi)^ ^EA;, see (j5.35j) . Hence, together we obtain, 
as i — > oo, 



Re(7f ^) + o(l) = Re ^7*'^ Y^EA*^ + o(l) 



(5.40) 



Thus, 7 7^ yields X^^^ E ^* = 1. For the bound on the Lipschitz constant in £2 of 
T restricted to see Rosier and Riischendorf [271 Lemma 1] and Fill and Kapur 
[nj: For choose (Z(i), W^(^)), . . . , (Z(^), T^(^)) as identically distributed 

vectors of optimal couplings of 11 and v and such that {Z^^\ W^-^^), . . . , {Z^^\ W^^^), 
(A^, . . . , A*j^) are independent. Then we have 



\r=l r=l / 



< E 



r) > 



r=l 



, r=l 



Ej^i^:i'^i(/i,i^)+o 

r=l 

YE\A;\'£l{^i,u). 



r=l 



Altogether we obtain that T has a unique fixed point C{X) in (t)- 
The fixed point property of C{X) implies 



1 Re (t^x) AlRe(j2 t^^lX^'^'^ 



(5.41) 



where {Al, . . .,A*^), X^^), . . . are independent and £(X('')) = C{X) for r = 

I, . . . ,K. We may assume, e.g. by taking optimal couplings, that \\Ar^ — A*\\2 —>■ 
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as t — > oo. We choose X^:^^ as optimal couplings to Re(t''^X('')) (with the right 
distribution, i.e. the distribution of Xt) for t > and r = 1, . . . , K . Clearly, we may 

(r) 

assume that, as required, X^ , r = 1, . . . , K , are independent of each other and of 
We denote, for t > 0, 



A{t) := l2 



Xu-Re[t^X 



Using and (|5.41|) we obtain, for t > tq, 



±r 1 



vr=l 



< 



< 



+ 



r=l 



+ 



+ 



bit) 



r=l ^ 



1 



X^*") ) Re ( t^AlX^''^ 



(5.42) 



By (|5.35() and ()5.33|) the second and third of the three latter summands tend to 
zero as t ^ oo. We abbreviate 



Hence, (|5.42|) implies 



X^ 2-1 Re 1 



1 



2\ 1/2 



(5.43) 



K 



K 



\ 



1/2 



V 



r=l 



r,s=l 
r^s 



+ 0(1). 



(5.44) 



/ 



By the definition of A(t) and the fact that (yX^\^e[t^'^ X^^^)) are optimal couplings 
for alH > and r = 1, . . . , we obtain 

2cr 



E(VF, 



(t)x2 



E 



(*)' 



t 



A2(tW) 



(5.45) 



From (|5..S2j) we obtain 



EXt = — Re(7t^) + R{t), t > 0, 
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with R{t) — > as t — > oo. Since KX^^^ = 7 and by the independence conditions we 



obtain EW^^^*^ = E[{T^'' /ty R{T)'')] and, for r / s, 



(^)^ 



i?(rW)i?(TW) 



Sphtting the latter integral into the events {Tr^^ < ti or t]*^ < ti} and {T^*^ > 
ti and t]*^ > ti} for some ti > we obtain, for every ti > 0, 



E[VFWm^W] < f ||i?||^ + sup 

V * / u>ti 



where ||-R||oo := sup^ < 00. From this we obtain first, letting t — > 00, 

limsupj^oo E[Wr^*-'iyi*''] < supu>j^ R^{u), and then, letting ti ^ 00, 



^ as t ^ 00. 
Now, (|5.44j) . (|5.45|) . and (|5.46j) imply, for t > tq, 

.(t) ^ 2a 



(5.46) 



A(t) < E 



E 



a2 ^ 



1/2 



+ +i?2(t), 



(5.47) 



with i?i(t), i?2(t) ^ as t ^ 00. 

We first show that ||A||oo < 00. Define A*(f) := supQ<„<j A(n). By the assump- 
tions supq<„<j E lYup < 00 and supo<tt<t |/(n)| < 00, together with = and 
f{u) = for u < 1, we have A*(t) < 00 for all t > 0. Let ti > tq be such that 
\Ri{t)\<l and |i?2(t)| < 1 for t>ti. Then with (I5l71l we obtain, for t>ti, 



A{t) < E 



1/2 

+ 11 +1. 



By (|5.33|) . (|5.35|) and (|5.36|) there exists a ^2 > such that for all t > t2 we have 
EE^i(^r*V*)^'' < ^ < 1. Thus, for all t > t2 we obtain, with ^/oT^ < ^/a + Vb 
for a, 6 > 0, 



and thus 



A(t) < yeA*(t) + 2, 



A*(t) < y^A*{t) + 2 + A*{t2) 



which implies || A||oo < (2 + A*(t2)) /(I - \/?) < 00. 

In a second step we show that A{t) — > as t — > cxd. For this we assume that 
L := limsupj^f^ A(t) > 0. Let e > 0. There exists a ts > t2 such that for all t > ta 
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we have A(t) < L + e. Then (|5.47j) imphes 



< E 



1/2 



r=l 



+ Ri{t)\ +R2{t) 



^ (j^^I^jY ml + ^{L + e)' + R,{t)^ +R2{t). 
Hence, t ^ oo imphes 

L<y^{L + e), 

which if L > is a contradiction if we choose e smah enough. Consequently, we 
have L = yielding the assertion. □ 

6. Proof of Theorem 11.31 

In this section we prove Theorem 11.31 The statements on mean and variance of 
N{x) are proved in Section |21 It remains to identify the asymptotic distribution of 
N{x) Note that recurrence (|1.3|) for N{x) is covered by the general recurrence for 
Yt in (|5.1() by making the choices d = 1, K = b, tq = 1, Ar{t) = 1, Tr^'' = Vr-t and 
bt = 1 for all r = 1, . . . , X and t > tq. 

We consider the three cases (i) - (iii) appearing in Theorem 11.31 separatelv: 

Case (i): We have EN{x) = a~^x + o{^/x) and Var(A^(x)) ~ j3x with /3 > 0. We 
apply Corollarv 15.21 with the choices f{t) = a~^t and g{t) = (3t. The conditions 
(|5.2U|1 and l|5.21jl are satisfied. We have sup„<^E|y„|* < cxd for s = 3 by Lemma 
Em Condition is satisfied with = ^/Vr for r = 1, . . . , /C, condition (031) 

is trivially satisfied, and we have 1)5. 24() . Hence, Corollarv 15.21 applies and yields 

N{x) — a~^x d 



which is the assertion. 



AA(0,1) 



Case (ii): We have EN{x) = a~^x + 0{^/x) and Var(iV(x)) ~ /3j;lnx with /? > 0. 
We apply Corollarv 15.21 with the choices f{t) = a~^t and g{t) = (3tlnt. Now we 
have g{T^:^^)/g{t) = Vr + Vrln{Vr)/ Int, hence we obtain, since x x In x is bounded 
on [0,1], 
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with A* = yJVr for r = 1, . . . , iC. All conditions of Corollary 15.21 are satisfied as in 
case (i) and we obtain 

^ px m X 

Case (iii): Theorem 15.31 can be applied and yields 

N[x) — a ^x jjjjj;^2inx^ 



X 



ReA2 



as x ^ oo. 

Here, we give a simplified version of the proof of Theorem 15.31 for the special 
recurrence H1.3|) which yields also the rate of convergence stated in Theorem 11.31 

The restriction of T defined in (|1.7j) to J^y\^{y^ is Lipschitz in I2 with Lipschitz 
constant bounded by (lEX]r=i Vr^^'^^^^^^'^ \ cf. the first part of the proof of Theorem 
15.31 Recall A2 = o" + ir with real o", r. Hence, o" > 1/2 implies that T has a unique 
fixed point in M^{v). 

For Xt := N{t) - a'H we obtain with (fO|) 

h 



^*=E (6-1) 



r=l 



(r) 

where are independent distributional copies of Xt also independent of 

(Vi, . . . , Vfc). With the fixed point property of H we have 

fRe (Se'"''^^) = Re{t^^E) = Re ^^(y.t)^^^^^ , 

where (Vi,...,Vfe), H^^^ . . . , H*^*) are independent and C{E^^^) = C{E) for r = 
l,...,b. We choose xj:^^ as optimal couplings to Re(t''*2H*^'')) for t > and 
r = l,...,b and denote A(t) := ^2(^*5 Re(t'^^H)). Note that in the definition of 
Xt we did not rescale by , hence we have to show A{t) = 0{t'^). 
With TyJ*^ := XQ - Re((V;.t)^2H(^)) we obtain, for t > 1, 

Ait) = £2 <.i + 1, E R« ((Kt)'^sM) ) 

\r=l r=l / 



2>| 1/2 

+ 1 



r=l r,s=l 



1/2 

+ 1. 



Conditioning on (Vi, . . . ,Vb) yields E(TVr'*^)2 = E A2(y,,t). From EN{t) = a'^t + 
Re(7t^2) + o^^K) and EH = 7 we obtain EW^^ = 0{t^). Since VfJ*^ and wP are 
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independent for r / s conditionally on [Vi, . . . , V;,), it follows that 



Mt) < l^EA'^iVrt) + 0{t^'') \ +1, t>l. (6.2) 

Now, we show that A{t)/t'^ = 0(1). Note that this implies the assertion. We denote 

^*{t) := sup 

l<u<t 

Then, (|(i.2j) implies, that for appropriate R > 

( b >j 1/2 

and, with a/o + 6 < ^ + \/6 for a, 6 > and ^ = E X]r=i < ^ ^^^^ implies 

**(t)<^-^<oo. 

The assertion follows. 



7. Examples 

Example 7.1 (Random splitting of intervals). Sibuya and Itoh studied the 
case of random splitting of intervals, with uniformly distributed splitting points; 
this is the case 6 = 2 and V = ([/, 1 - U), with U ~ U(0, 1). 
We have 

(t){z)=^U^ +^{l-UY = 2 ! u'du = -^, Rez>-1, 

Jo I + z 

which is a rational function. The characteristic equation p.5() is 2/(1 + A) = 1, 
and has the single root A = 1. Thus Theorem 11.3( 1) applies and shows asymptotic 
normality, as stated by Dean and Majumdar 7 . Further, a = — (/>'(1) = 1/2, so 
Theorem l^.ir ii) yields E A^(x) = m{x) = 2x + 0{x^) for every 6 > 0. More precisely. 
Theorem 13.41 yields 

EiV(x) = m(x) = 2x - 1, X > 1, 
which also can be shown directly from H1.2|l or from ()3.11|) . 
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For the asymptotic variance, we obtain from Theorem l,'-{.4r ii) . since M = 1 and 
ao = —1, using symmetry, 

P = a-^ (e{U + 2C/ a (1 - [/) + 1 - [/) - 2 + 1^ 

-31 mwfTT A /I TT\\ 1 „.— 1 



+ a-" [2E{U A {1 - U)) j - a 

/■1/2 

= 20E([/ A (1 - U)) + 16 / n(ln(l - u) - In('u)) du - Q 

Jo 

= 81n2 - 5 w 0.545177. 
This can also be obtained from Theorem 13. If iii) : we have 

i;{z,w) = e(^{U' + (1 - Uy)^^" + (1 - U)"")^ - 0(z)0H 

^ +2B{z + l,w + l) ^ 



l + z + w ' ' ' {l + z){l + w) 
2 ^^T{z + l)T{w + l) 4 



l + z + w T{z + w + 2) (l + z)(l + u;) 

and thus 

4 

-0(1/2 + iu, 1/2 - iu) = 1 + r(3/2 + i'u)r(3/2 - iu) 



|3/2 + in|2 



1 + |l/2 + in|2 



coshvru |3/2 + iuP ' 

and, since 1 — <^(z) = {z — l)/{z + 1), 

n f TT , , ,2 4 \ |3/2 + i'u|2 , 

/3 = -/ 1 + — ^ 1/2 + iuP--- ^ ^dn, 

vr^.o^V coshvrii' ' ' |3/2 + inp/ |l/2 + i-uj^ ' 

which can be integrated (with some effort) to yield 8 In 2 — 5. 
Consequently, by Theorem II. 3| 

^(-^-^-^AA(0,81n2-5). 

Example 7.2 (m-ary splitting of intervals). We can generalize Example 17.11 bv 
splitting each interval into m parts, where m > 2 is fixed, using m — 1 independent, 
uniformly distributed cut points in each interval. This has been studied by Dean 
and Majumdar 

We have b = m, and Vi, . . . , Vm have the same distribution with density (m — 
1)(1 - x^-^, < < 1. Hence, 

(f){z) = mKVi = m{m — 1) / x'^ {1 — xY dx = m{m — l)B{z + l,m — 1) 

Jo 

T{z + l)m\ ml 
r(m + z) {z + 1) ■ ■ ■ {z + m — 1) 
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The characteristic equation (j){z) = 1 becomes T{z + m)/T{z + 1) = ml, or 

{z + !)■■■ {z + m- 1) = ml. (7.1) 

The same equation appears in the analysis of m-ary search trees. It is shown by 
Mahmoud and Pittel [121 and Fill and Kapur [12] that if m < 26, then ReA2 < 1/2, 
and thus (i) applies, but if m > 27, then ReA2 > 1/2, see also, e.g., Chauvin and 
Pouyanne [31 and Chern and Hwang [Hj. Theorem K14I yields an exact formula for 
EiV(x) (although it is hardly useful except when m is small). It further leads to a 
formula for the asymptotic variance, provided m < 26. 

We have, with ip{z) := r'{z)/T{z) and Hz := ip{z + 1) — (for integer z, these 
are the harmonic numbers) 

a = -0'(1) = ij{m + 1) - ^(2) =H^-1. 

Example 7.3 (Random splitting of multidimensional intervals). Another general- 
ization is to consider d-dimensional intervals, where an interval is split into 2'^ sub in- 
tervals by d hyperplanes orthogonal to the coordinate axis and passing through a 
random, uniformly distributed point. This too has been studied by Dean and Ma- 
jumdar [3|. 

We have b = 2'^. Vi, . . . ,Vb have the same distribution, Vj = Ui - ■ ■ Ud, where 
Uk ~ U{0, 1) are i.i.d. Hence, 

(l){z) = 2'^EV{ = 2'^{EU!Y = (r^)''- 

Again, (f> is rational. The characteristic equation may be written ((1 + A)/2) = 1, 
with the roots 

A = {26^^'''^'^ - I : < k < d - I}. 

Thus CJ2 := ReA2 = 2cos^ — 1, and the condition ReA2 < 1/2 is equivalent to 
cos(27r/(i) < 3/4, which holds for d < 8, while Re A2 > 1/2 for d >9. This justifies 
the claims in Dean and Majumdar j7j. 

The same characteristic equation, and the same phase transition, appears for 
quad trees, see Chern, Fuchs and Hwang 1^. 

We further observe that a = —(f)' (1) = d/2. 

The random trees in these three examples have also been studied by ^^1; QH] 
and ^nj) where the properties of a randomly selected branch are investigated. This 
problem is quite different, and there is no phase transition. See also |17j . 

Example 7.4 (Random splitting of simplices). Consider d-dimensional simplices, 
where an simplex is split into d + \ new simplices by choosing a random point X 
in the interior and connecting it to the vertices of the original simplex; each new 
simplex has as vertices X and d of the original d+1 vertices. 
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It is easily seen that this is equivalent to d+l-ary splitting as in Example 17. 2( see 
[SI Lemma 3], so we have the same results as there, with m = d + 1. In particular, 
N{x) is asymptotically normal if < 25. 

Example 7.5 (Non-uniform splitting of intervals). Returning to binary splitting 
of intervals, we can generalize Example 17.11 by taking another distribution for the 
cut points; we thus have b = 2 and V = {V, 1 — V), where V has any distribution 
on (0, 1). An interesting case is when V has a beta distribution V ~ B{a,a') with 
a, a' > 0; then 

Jo B{a,a') T{z + a + a') T{a) 

E(l — Vy is obtained by interchanging a and a' . In particular, if a and a' are 
integers, then cp is rational. 
We consider two special cases. 

(i) The symmetric case with a' = a, V ^ B{a,a). Then 

S(z) - 2 ^(^ + ") ^(^"^ - r(z + a) r(l + 2a) 
^^^^ ~ r{z + 2a) r(a) ~ r(z + 2a) r(l + a) ' 

We have a = —(p'{l) = H2a — Ha, with as in Example 17.21 Numerical 
solution of the characteristic equation seems to show that ReA2 < 1/2 if 
and only if a < cq, where ~ 59.6. 

(ii) The case a' = \, V ^ B(a, 1). Then 

Tjz + a) r(a + l) , T{z + 1) a , r(z + 1) 

(p{z) = —, r — ——^ ^— -1 a+1 = \-—, ri a+1 . 

^ r(z + a + l) r(a) r(z + a + l) ^ ^ z + a T{z + a + l) ^ ' 

One finds a = Ha/ {a+1). The characteristic equation </>(A) = 1 is equivalent 

to r(a + l)r(A + l)/r(A + a + l) = A/(A + a) or 

r(a + A) 

When a = m is an integer, this is the same as ()7.1() . so Re A2 < 1/2 for integer 
a if and only if a < 26. In general, numerical solution of the characteristic 
equation seems to show that ReA2 < 1/2 if and only if a < oq, where 
ao ^ 26.9. 

8. NON-EXAMPLES 

In this section, we give a few examples where our theorems are not valid. 

Example 8.1 (Lattice). In the lattice case, there exists R > 1 such that every 
Vj G {R-^ : /c > 1} U {0} a.s. In this case, (j) is periodic with period 27ri/lni2; in 
particular, the characteristic equation H1.5() has infinitely many roots 1 + 27rin/lni? 
on {A : Re A = 1}, and thus Condition B((5) fails. Indeed, it is obvious from (|1.2j) 
that N{x) = N{R'^) when R"^ < x < so EiV(x)/x oscillates and does not 
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converge as x ^ oo. The natural approach is to consider only x G {i?™ : m > 0}. 
It is then straightforward to prove an analogue of Theorem 13.11 using the lattice 
versions of the renewal theory theorems that were used in Section |21 An analogue 
of Theorem II .31 then follows by the usual (discrete) contraction method, as in ^4'. 
We leave the details to the reader. 

Example 8.2 (Deterministic). lfV = (Vi, . . . , Vb) is deterministic, then so is N{x), 
and it is meaningless to ask for an asymptotic distribution. However, it makes sense 
to study the asymptotics of N{x) = m{x). (Clearly, (T^(x) = 0.) 

If V is non-lattice, then N[x)/x ^ a by Theorem 13.11 and Remark 13.21 If V is 
lattice, we consider, as in Example 18.11 only x = R"^, m >1. 

We may assume that V^- > for each j. By Dirichlet's theorem jl4l Theorem 
201], for every e > 0, there exist arbitrarily large t such that \Vj^ — 1| < e for 
j = 1, . . . ,b; thus limsupf_,oo \<j){l + it)\ = 1. Hence Condition B(l) does not hold, 
and therefore, by Lemma l2.H Condition B(J) does not hold for any 5 < 1. 

More precisely, if | — 1| < e for j = 1, . . . , 6, let = 1 + it. Then \(j){zQ) — 1| < e 
and 

\(p'izo) + a 

Since further |(/)"(z)| < | InV^ p for Rez > 0, it follows easily that if e is small 
enough, then (j)(z) — 1 has a zero in the disc B := {z : \z — zo\ < 2e/a}. (Use the 
Newton-Raphson method, or Rouche's theorem and a comparison with the linear 
function 4>{zo) + {z — ZQ)(j)'{zo).) It follows that there exists a sequence A„ G A with 
Re An ^ 1 and ImA^ — > +oo. 

We give some concrete examples: 

V = (1/2, 1/2) is lattice with R = 2 and iV(2") = 2". 

V = {t^^jT^"^) where r = (1 + \/5)/2 (the golden ratio) is lattice with R = t 
and N(t'"') = -Fn+a — 1, n > 0, as is easily proven by induction. (F„ denotes the 
Fibonacci numbers.) Thus, A^(r") ~ 5~"'^/^r"+^. 

V = (1/3,2/3) is non-lattice and thus N{x) ~ a~^x, where a = ^ln3 + 
|ln(3/2) = ln3 - |ln2. 

9. Some related models 

The basic model may be varied in various ways. We mention here some variations 
that we find interesting. We do not consider these versions in the present paper; we 
leave the possibility of extensions of our results as an open problem, hoping that 
these remarks will be an inspiration for future research. 

Remark 9.1. By our assumptions, the label of a node equals the sum of the labels 
of its children. Another version would be to allow a (possibly random) loss at each 
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node. One important case is Renyi's parking problem ,26,, where a node with label 
X is interpreted as an interval of length x on a street, where cars of length 1 park at 
random. Each car splits an interval of length x > 1 into two free intervals with the 
lengths U{x — 1) and (1 — — 1), where U ~ U(0, 1). An obvious generalization 
is to split {x — 1) using an arbitrary random vector (Vi, . . . , V;,). (The one-sided 
version, where we study only one branch of the tree, is studied in JHl; |17j.) 

Remark 9.2. Krapivsky, Ben-Naim and Grosse j211 l2Uj have studied a fragmen- 
tation process where fragmentation stops stochastically, with a probability p{x) of 
further fragmentation that in general depends on the size x of the fragment. Our 
process is the case p{x) = l{2;>i}. Another interesting case is p{x) = 1 — e~^, see 
Remark 19.31 below. 

Remark 9.3. Our model is a continuous version of the split trees studied by De- 
vroye ^8,, where the labels are integers (interpreted as numbers of balls to be dis- 
tributed in the corresponding subtree) and each label n is, except at the leaves, 
randomly split according to a certain procedure into h integers summing to n — sq; 
here sq is a small positive integer (for example 1) that represents the number of 
balls stored at the node. Typical examples are binary search trees, m-ary search 
trees and quadtrees. We can regard the continuous model as an approximation of 
the discrete, or conversely, and it is easy to guess that many properties will have 
similar asymptotics for the two models. This has been observed in several examples 
by various authors, see |7j and For example, the results for Example 17. 21 parallel 
those found for m-ary search trees by |^, [5], ^U, |3I and others. Similarly, the 
results in Example 17.31 parallel those found for quadtrees by • 

We study only the continuous version in this paper. It would be very interesting 
to be able to rigorously transfer results from the continuous to the discrete version 
(or conversely); we will, however, not attempt this here. 

Note that for binary search trees, we have n random (uniformly distributed) 
points in an interval, split the interval by the first of these points, and continue 
recursively splitting each subinterval that contains at least one of the points. If we 
scale the initial interval to have length n, then the probability that a subinterval 
of length X contains at least one point is ~ 1 — . Thus it seems likely that 
the binary search tree is well approximated by a fragmentation tree, with V as 
in Example 17.11 with a fragmentation probability 1 — as in Remark 19.21 The 
same goes for random quadtrees and simplex trees corresponding to Examples 17.31 
and El 
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